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Abstract 

Generation mechanism of energy gaps between conductance and valence bands 
is at the centre of the study of graphene material. Recently Chamon, Jackiw, et 
al. proposed a mechanism of using a Kekule distortion background field (p and its 
induced gauge potential Ai to generate energy gaps. In this paper various vortex 
structures inhering in this model are studied. Regarding (/? as a generic background 
field rather than a fixed Nielson-Oleson type distribution, we have found two new 
types of vortices on the graphene surface — the velocity field vortices and the 
monopole-motion induced vortices — from the inner structure of the potential A^. 
These vortex structures naturally arise from the motion of the Dirac fermions instead 
of from the background distortion field. 
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1 Introduction 



Graphene is a lately realized single-atomically thin 2- dimensional sheet of carbon. Its non- 
Bravais honeycomb lattice arrangement of atoms gives rise to a new many-body problem 
that single-particle quantum mechanics is relativistic, but interactions are non-relativistic 
instantaneous Coulomb interactions. Ever since the breakthrough of its manufacturing 
technique in 2004 [1] graphene has attracted enormous interest because of the funda- 
mental new physics exhibited and the highly unusual electronic, mechanical and optical 
properties. These properties lead to numerous potential application and make graphene 
a new material of growing technological importance [2] . 

Graphene has a remarkable electronic property of ballistic transport at room tem- 
perature under electrical and chemical doping. This makes graphene a candidate to be 
the next dominant semiconductor material in electronics after silicon. However, a major 
obstacle for this application is that the energy spectrum of monolayer graphene remains 
metallic at the Dirac points because of the relativistic linear dispersion relation, i.e., the 
vanishing gap between the conductance and valence bands in the energy spectrum. Many 
methods have been developed for generating a non-zero gap, such as spatial confinement 
and lateral super-lattice potentials. In this paper we mainly consider the generation mech- 
anism proposed by Ghamon, et al. [21 HJ O [6l [TJ El |9] and aim to find new topological 
vortex structures originating from the singular configurations of the Dirac fermion field. 

The Hamiltonian for the honeycomb lattice of monolayer graphene is given by 



3 

Ho = -tJ2Yl ^ + + bH^ + Si) a (x)] , (1) 



xeA i=l 



where a and b are respectively the fermion operators acting on the sublattice A and B, 
and t the uniform hopping strength. In the momentum space Eq.([T]) is written as 

ifo = [r (k) at (k) b (k) + r* (k) b^ (k) a (k)] , (2) 

k 

{a (k) I _., I a (x) 1 3 ., 

> = 2Jx^ «Kx y y g^j^j TCk) = —t22^^^e '. The single-particle 

6(k) J [ 6(x) J 

energy spectrum contains two Dirac points K-t = ± (^^^;0j at the zone boundaries, i 
being the length of lattice, satisfying t(K±) = 0. In the neighborhood of K-t, k is written 
as k = K-t + p and Hq is linearized as 

Ho = $^[r+(p)aV(p)6+(p) + r;(p)6V(p)a+(p) 
p 

+r_ (p) al (p) b_ (p) + T*_ (p) bl (p) a_ (p)! , (3) 



where a± (p) = a (K± + p), 6± (p) = b (K± + p) and r± (p) = r (K± + p) = ±vf {px ± ipy), 
Vp being the Fermi velocity, fi? = 1. Hq has the following matrix form: 

Ho = 1 = 1,2. (4) 

Here \E' (x) is a Dirac 4-spinor, \1' (x) = ^ (x) , (x) , tp": (x) , (x) j , with 



(p) J I v^i (xo J 



. The cti, i = 1, 2, is given by = 7471, accompa- 



nied by the definition of the Dirac matrices in quantum electrodynamics: 7^ 



-(Ik 

(Tfc 



/c = 1, 2, 3; 74 = /3 = ( ^ J ; 75 = "^74717273 = ^ ^ 1 ; "^fc Pauli matrices. The 

7-matrices are Clifford algebraic 1-vectors, satisfying 7^76 + 7f,7a = 26abl, a,b = 1, 2, 3, 4 
[TU] . From (jl]) one obtains the equation of motion for : 

7a9a^ = 0, a = 1,2,4, (5) 

which is recognized to be the massless Dirac equation in (2 + l)-dimensions. The energy 
dispersion of (jSD is a linear relation: e (p) = ± |p|, which has vanishing gaps between the 
conductance and valence bands in the energy spectrum. 

For the purpose of generating energy gaps, Chamon, Hou and Mudry [ij |3] introduced 
to the monolayer graphene a complex scalar background field A (r). The A (r) describes 
a Kekule distortion of the lattice and leads to a small variation 6tr,i over the unit hopping 
strength t : 

5tr,i = [A(r)e^^+-^'e*^-'" + A*(r)e~*^+-^'e-*^""] , G = K+ - K_, (6) 

Stj- i yielding chiral mixing between K+ and K_. With ([6]) the Hamiltonian ([1]) is modified 

as 

3 

H = -J2J2^^ + '^^r,) [at (r) b (r + s,) + 6+ (r + s,) a (r)] , (7) 

reA i=l 

and the matrix form (jlj) becomes 

H = $t/c^ = $t [«.^^ + ^^^^ _ ^, ^ = 1, 2, (8) 

where the kernel /C is 



( -2id-, A(r) \ 

-2i9j A(r) 

A*(r) 2id^ 

\ A*(r) 2i(9j / 



(9) 



In ([H]) one denotes A (r) = gip, where g is a. coupling strength describing the dimension of 
A (r), and ip = v^Re+^V^im a complex scalar, y^^e and ipi^ being the real and imaginary parts 
respectively, can alternatively be expressed via its modulus and phase as ip = \ip\e^^. 
The energy level equation derived from (|8]) reads 

[aiPi + gl3 {ifRe - hmm)] ^ = E^. (10) 
It is pointed out IH [5] that x can be removed by a gauge transformation of ip : 

= ^^ + 2uj, ^ ^ = 6^"^^^^. (11) 

Thus dH]) yields {'jada + iglfl)"^' = 0. This gives rise to a mass in the dispersion relation, 
e (p) = i-^/p^ + I A|^, and hence generates a desired energy gap between the conductance 
and valence bands. 

For the purpose of keeping the kinetic portion of ffTOj) invariant against the local 
gauge transformation ( ITTj) . Jackiw and Pi ^ extended (fTOj) by introducing a U (1) gauge 
potential Ai : 

[a,{p,-^5A)+9HfRe-i^5V>im)]'^ = E<i/, 2 = 1,2, (12) 

where Ai transforms as Ai — Ai + diU under f lTTj) . (In f fT2|) the chemical potential is not 
considered |6].) This new model f|T2|) is an equation about the Dirac fermion wavefunction 
\& and the background field ip and its induced gauge potential Aj. In [5j \l/ is solved out 
by fixing Ai and (p within the Nielsen-Olesen (N-0) vortex configuration. Very recent 
developments based on the above model include the study of the superconducting vortices 
by P. Ghaemi, S. Ryu and D.H. Lee [7] and P. Ghaemi and F. Wilczek [8], and the study 
of the fractionally charged vortices by J.K. Pachos, M. Stone and K. Temme [S]. 

In this paper we will present another analysis for ip and Ai. Firstly, our starting point 
is that the background field ip permits more generic distributions on the graphene surface 
rather than a fixed N-0 type, since ip is an externally given distortion field. Correspond- 
ingly, we regard Ai and as two unknown quantities that can be determined by the 
governing equation (fT2|) together with supplementary restraints, where the restraints are 
naturally provided by the subsequent study of topological vortices. Secondly, two other 
types of topological vortex structures can be derived from the physical equation (fT2|l — 
the velocity field vortices and the monopole-motion induced vortices. These vortices arise 
from the singular configurations of the spinor field \1/ instead of from the background ip. 
Thirdly, moreover, it will be shown that the potential Ai contains a term which describes 
the direct coupling of ip and \1/ fields and can be expressed formally by an 5*0 (4) gauge 
potential. 

This paper is arranged as follows. In Sect. 2 the inner structure of the gauge potential 
Ai in terms of the Dirac fermion spinor wavefunction \1/ is studied. In Sect. 3 two types of 



topological vortex structures in the graphene system are studied. In Sect. 4. the paper is 
summarized. 



2 Chiral Gauge Theory for Graphene 

In following it is shown that introducing the Ai in ( !T2|) implies formally introducing one 
SO (4) potential and two U (1) potentials. 

Consider the zero-energy level of the model (JT2i) : 



7i {di - ij^Ai) ^ + ^ {iipKc + ^ = 0, i = l,2. (13) 

Left-multiplying (fT3|) with \E'^7j75 there is 

^Wi757i<9i^' + iAi^H^-fj-fi^ - i^¥;Re^^"^757j^ + gcpi^^-fj^ = 0. (14) 

f[T^ minus its complex conjugate gives a decomposition expression for in terms of the 
spinor wavefunction "if [HI [T2] : 



A = Af^ + Af + Af, (15) 



where 



\]>t\[/ vl>t\[/ 



(16) 



being the SO{A) generator, Uj = | [7j,7j] • Eq.(lT5l) is a sufficient but not necessary 
condition for the original equation (!T3|) . In this section and the next the three terms 
Af \ Af^ and ^4^^^-' will be investigated respectively for their topological and gauge field 
theoretical essence. 

A-^'' contains the coupling between \1/ and the background distortion field (p. The 
discussion is twofold: 

Firstly, Af '' can be formally expressed by an SO (4) gauge potential. Indeed, singling 
out the A[^^ part in f[T51) leads to 

jidi^ - tja^Af^^ + tgvKe'^ + 9^1^15"^ + Z757i (4'^ + Af) ^ = 0, t = l,2, (19) 

which can be compared to the following equation 

'yA^ -laUJa^ + ^9^Rc^ + 9^lml5^ + ^75li (^f ^ + ) ^ = 0, o = 1,2,3,4. (20) 



Here Ua = \(jJabchc is an 5*0 (4) potential, ujabc being anti-symmetric for the h c 
indices. Our aim is to express f|T9l) by fl20|) . hence we write Uabc in three parts: 
^abc = ^abc + ^'abc + ^abc^ whcrc cu^j,^ is fully anti-symmetric for a b c, oo^^^ symmetric 
for a b, and u^^^ symmetric for a c: u^^^ = ^iuabc + ^bca + ^cab), ^llc = |(waf,c + Wfcac), 
^abc = li^abc + (^cba)- Defining Uh = ^Uaba, = (^abcd^td^ One obtains after simple Clif- 
ford algebra: ujabd^^lblc = hal^^^a, ^lldalblc = -|c^c7c, ^Ibclalblc = -|c<;67fe- Then 
( 120|) becomes 

jidi^ - ^7, (275W, - a;,) + igy^Ke"^ + g^iml^"^ + il^li [^^? + Af^ * = 0. (21) 

Thus comparing ( l2Tl) with ( !T9l) one arrives at a conclusion that if the evaluation 

uj^ = 0; bJi = 4Af \ « = 1, 2; cjg = CU4 = (22) 

is taken, the 2-component potential a!>^^ ^ i = 1,2, can be formally expressed by the 4- 
component SO (4) potential Ua, a = 1,2,3,4, via the expression ( l20l) . This is the gauge 
field theoretical essence of the potential A^^^ It is addressed in addition that since Af^ 
is real and 7^ anti-Hermitian, the 5*0(4) potential Ua is required to be Hermitian, which 
is different from ordinary anti-Hermitian gauge potentials in physics. For the discussion 
of Hermitian gauge potentials see [I3] . 

Secondly, Ai is the induced gauge potential of the background distortion (p, and A^^^ 
the only term of Ai which explicitly contains (f. Our second observation for Af'^ is that 
(f enters the potential A^ by coupling to the Clifford algebraic odd vectors [TU]. Go back 
to ([13]) and write the Af^ part explicitly 

+ ^757i (^f + ) ^ + ^^7<^Re^ + mm75^ = 0, ^ = 1,2. (23) 

Keep in mind that 7^ is a Clifford algebraic 1-vector and 7574 a 3-vector (or pseudo 1- 
vector). Then introducing the projection operator P = Eg. (1231) reads 

-fidi'^+g (i^nePr.ir + <^imA»75) ^+«757i {Af^ + ^) ^+ic/(^Re^ + C/(^I„75^ = 0. (24) 

Here -Ryg-y^ = Tr (^Pj^ji^ 7571 and P-^- = Tr (^Pli^ 1% are the Clifford vectorial components 
of P respectively along the 3-vector 757^ and the 1-vector 7j. Hence it is concluded from 
(!24|) that if enters the potential A^ by coupling to the Clifford odd vectors. 

Moreover, it can be checked that A^-^^ and Af^ are invariant under the gauge transfor- 
mations (ITT!) of Lf, whereas Af''^ undergoes a f/(l)-type gauge potential transformation: 
a!? ^ 4" = ^\ A^ ^' = Af^ A^' ^ A^" = A^' + d.u. 



3 Topological Vortex Structures 

In this section the Af'^ and A^f^ terms of ( IT5|) will be investigated. It will be shown that 
they lead to different types of topological vortices on the graphene surface. 

(2) 

Vortices Arising from Al 

With respect to the chiral 75 matrix one considers the self-dual and anti-self-dual part 
of separately: = + with vE^^ = i (1 ± 75) ^ satisfying 75 vE^^ = ±^1^^. The 

explicit representation of \E'± is: \I'+ = ^ $^ j , = ^ , where $+ (x) 
and $_ (x) are two Pauli 2-spinors and have their respective SU (2) sub-group spaces 
(hereafter marked as SU (2)_,_ subspaces). Thus (fT7|) reads 

(2) ^ ^^^2 ^■'^+^+ - sin2 

2i$^+$+ 2z$l$_ 

where cos^ 6^ = and sin^ 9 = ^^^^ , 6* (x) called the duality rotation [TT] . 

As mentioned in the introductory section, in the present paper the gauge potential Ai 
and the spinor wave function are regarded as two unknowns determined by the governing 
equation f|T2|) . with the background field cp being a generic externally given distortion for 
the graphene system. f[T^ provides 8 restraints but Ai and \E' have totally 10 unknown 
components. Hence for determining the unknowns one needs two extra supplementary 
restraints. In this regard one employs the following normalization conditions of \E' to play 
the role of the required supplementary restraints: 

sin^ 9 = cos^ 9 = ^, = 1. (26) 

Eq.( l26|) means that there is no priority between the self- and anti-self-dual subspaces, 
hence (|26l) is a natural choice. Thus 

= lll ' = - H^^"^^ - • (27) 

(2) (2) 

The bilinear form A]^ (x) [resp. A-_ (x)] is a f/ (1) gauge potential in the SU (2)^ [resp. 
SU (2)_] subspace, and has the physical meaning of the velocity field of spinning quantum 
fluid, up to a mass constant. 

Topologically singular property of the configuration of $-1- is studied by considering 
the first Chern class, a topological characteristic class on the 2-dimensional graphene 

('2) f2'l 1 1 ('2'l 

surface [H]. The first Chern class constructed from Ai^ is C^J. = ■^^f-j^dx'^ Adx^ , where 
fijl- is the gauge field strength defined as f^Jl. = diA^^_l — OjAI'^. The physical meaning 

(2) 

of C{_(_ is the vorticity of the velocity field distributed on the graphene surface. It is 



known in hydrodynamics that non- vanishing vorticity possesses vortices; this fact can be 
geometrically shown here for Cl^. Indeed, a spinning quantum fluid such as the Helium-3 



superfluid has the so-called Mermin-Ho relation [15] 

d.Afl - djAf2 = leabcnld,nidX±d^x, (2^ 



where n^., a = 1, 2, 3, denote the spin unit vectors respectively in the SU {2)_^ subspaces, 
(x) = ^r^^- From the geometric point of view, the RHS of f l28p is the pullback of 
a surface element on the hypersphere S'^, where 5*^ is formed by in the group space 
of SU {2)_^. The tangentially intersects the graphene surface at the studied point 
X = (a;^,x^), hence is locally perpendicular to the graphene surface at x. Meanwhile, 
it is known that the S'^ surface element can be geometrically given by a f/ (l)-type gauge 
potential w± as 

eabcn'^dn'^^dn'^ = dw±. (29) 

Here w± is defined as w± = ei±-de2±, where ei+ and 6*2+ [resp. ei_ and e2_] are a pair of 2- 
dimensional unit vectors normal to 77,^^ [resp. n°i]onS'^ [resp. S"^]. Namely, (ei+, 6*2+, ri+) 
[resp. (ei_, e2_, n_)] forms an orthogonal frame: e±i ■ e±2 = e±\ ■ n± = e±2 ■ n± = 0, 
e±i ■ e±i = e±2 ■ e±2 = n± ■ n± = 1. The w± has the physical meaning of the so-called 
Wu-Yang potential ^j6j. Comparing ( |28|) and (p9|l one knows that locally 

A± = 'Wi± = ei± ■ die2± (30) 

up to a removable phase angle. Hence in following one can use fl30p to investigate the 
topological structures arising from A]^. 

For convenience we use one unique 2-component vector to replace the pair (ei±,e2±). 
This vector, denoted as ^± = {^]_,^'^), is required to reside in the plane spanned by ei± 
and e2± and satisfy ef^ = = e^-^^, with = ^^^^ and A, B = 1,2. The 

zero points of ^± are the singular points of ei± and e2±- Then in terms of C,± one rewrites 
A^^2 = Cab iii^c^j and then obtains = ■^eijeABdij^^^dj-^^^d'^x. According to 
it can be proved that 

C[2 = 5' (el) D{^-^)d\ (31) 
where D{C,±/x) = \e^^ eAsdi^^dj^^ is the Jacobian determinant. Eq.f l3ip shows that non- 

(2) 

vanishing ClJ. occurs only at the zero-points of ^-i-. Hence in order to find the expected 
topological vortices on the graphene surface one should study the zero point equation 
of ^^{x) = 0, A = 1,2. The implicit function theory [TS] declares that under the 
regular condition D{^±/x) ^ the general solutions of the zero point equations are a 
finite number of 2-dimensional isolated points: 



(32) 



where L± denotes the number of the isolated points. Thus these points are the predicted 
topological singularities in the configuration of the $-|- field on the graphene surface. In 
the scenario of hydrodynamics they are called the velocity field vortices, the topological 
vortices arising from non- vanishing vorticity. 

The (5^ (^) ^^"^ expanded onto these vortex points as 

5^ {i±j = J2k±=i ^k±S^ {x — Xk^) , where Wk^ is the winding number of the fc-t-th zero- 
point, playing the role of the topological charge of that zero-point. Therefore the first 
Chern number given by the Chern class C^^ is 

cfl = I Cf2 =J:W,,. (33) 

k± = l 

(2) 

has the physical meaning of quantized vorticity. 

(3) 

Vortices Arising from A\ 

Another kind of vortex structures originate from the third term of Eq. ( fT5|) , Af^ = dj [\E'"''75/jj\E'] . 
Defining the dual tensor of lij as = 75/ij = ^eijkihi {hj = 1,2; k,l = 3, 4), one has 

Af = [^^ . (34) 

The RHS of f p4|) can be expressed by a Maxwell-type U (1) field tensor. Indeed, according 
to [12] a Maxwell-type U {1) electromagnetic field strength F^y has a Dirac spinor repre- 
sentation F^y = \l/"'"J^i,\l/, /i, z/ = 1,2,3,4, provided that \1/ is non-singular. Therefore, for 
( 134|) one can extend the studied base space from the 2-dimensional to the 4-dimensional, 
and consider only the weak solutions for Here the weak solutions refer to the case that 
the \l/ field has only a countable number of isolated singular points and is well-defined 

almost everywhere except at the singular points. Then the RHS of f l34p can be expressed 

1 , 



by a dual electromagnetic field strength as a[^^ = dj *Fji, where *Fij = \eijkiFki is the 



dual tensor of Fij. 

The second Maxwell equation reads dj *Fji = —4tt *Ji, where *Ji is the current of 
monopoles. When * Jj = 0, namely there is no existence of monopole-type excitations in 
the system, there is dj *Fji = 0, which corresponds to the Bianchi identity. In this case, 
Ai ' =0; When *Ji^O,A\' = —Att *Ji. One needs to consider the vortex structures 
arising from the motion of monopoles on the 2-dimensional graphene surface. 

There are various models for magnetic monopoles in literature [Til [T^ . depending 
on different choices for the symmetry of the monopole field or the gauge potential. In 
following we simply consider the complex scalar field model for the monopoles, and let 
the monopole current and density take the following form 

*J^ = ^r^{r^^^P-^iri^), 1 = 1,2- p = V>, (35) 



where ip is a complex scalar, ip = (f)^ + i(f)^. ( 135|) leads to *Ji = e""^ (^) ^* (w) ^^^^ 

= a, /3 = 1, 2. The first Chern class induced by Af^ is ^ = le^'^e'^di (^g) 6»j 

Since only the case of weak solutions is taken for ,C\ can be expressed in the 5- function 
form as c'f'^ = 5^ (^(p^ D (^), similarly to the last subsection. Apparently non- vanishing 
C\ occurs only at the zero-points of 0, so the zero point equation = 0, a = 1,2, 

should be studied. Similarly, under the regular condition D (0/x) 7^ the general solu- 
tions of the zero-point equations are a finite number of 2-dimensional isolated points: 

x' = x\, £ = l,2,...,Ar, (36) 

where N denotes the number of the isolated points. These points are topological singular- 
ities in the configuration of the field on the graphene surface; they are vortex structures 
due to the motion of monopoles. Furthermore, the 5^ can be expanded onto the vor- 
tex points as 5^ = J2f=i Wi6^ (x — Xi), where Wi is the winding number of the £-th 

zero-point, playing the role of the topological charge of that zero-point. Correspondingly 

(3) 

the first Chern number given by the Chern class is 

.(3) 



J Cf ) = J2 (37) 



4 Conclusion 



The massive Dirac equation ( 1T5|) proposed in Ref. [5] is an important mechanism for gener- 
ating energy gaps between conductance and valence bands of graphene material. We start 
from this equation and proposed that the distortion field (f is permitted to take a generic 
distribution rather than a fixed Nielson-Oleson (N-0) configuration on the graphene sur- 
face. The (y9-induced gauge potential Ai and the Dirac fermion wavefunction are de- 
termined by f|T3|) and the normalization condition fl26|) . In this paper we focus on the 



inner spinor structure of Ai and study its three terms ^4^^^'', Af'^ and Af\ Emphasis is 
on two new types of topological vortices arising from A] and A] — the velocity field 
vortices and the monopole-motion induced vortices. These vortices are different from the 
N-0 vortices, as they both originate from the singularities of the wavefunction \l/ instead 
of from the background field (p. Moreover it is shown that the A[^^ term can be formally 
expressed via an SO (4) gauge potential Ua- 
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